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We consider a thin elastic sheet adhering to a stiff substrate by means of the sur- 
face tension of a thin liquid layer. Debonding is initiated by imposing a vertical 
displacement at the centre of the sheet and leads to the formation of a delaminated 
region, or 'blister'. This experiment reveals that the perimeter of the blister takes 
one of three different forms depending on the vertical displacement imposed. As 
this displacement is increased, we observe first circular, then undulating and finally 
triangular blisters. We obtain theoretical predictions for the observed features of 
each of these three families of blisters. The theory is built upon the Foppl-von 
Karman equations for thin clastic plates and accounts for the surface energy of 
the liquid. We find good quantitative agreement between our theoretical predic- 
tions and experimental results, demonstrating that all three families are governed 
by different balances between elastic and capillary forces. Our results may bear on 
micrometric tapered devices and other systems where elastic and adhesive forces 
are in competition. 

Keywords: Adhesion; elastic plate theory; elastocapillarity; developable 

surfaces; buckling 

1. Introduction 

Adhesion is ubiquitous in a range of industrial applications and biological situations. 
On the one hand, highly optimized adhesives are used to assemble components or 
to repair broken objects. In electronics and in coating applications, an elastic film 
often adheres to a substrate through molecular forces (see Xia et ai, 1999, for 
example). On the other hand, geckos and insects use a combination of Van der 
Waals and capillary forces to adhere to substrates and walk upside down (Autumn 
et ai, 2002; Huber et ai, 2005; Gorb, 2005). 

While industrial applications generally seek permanent bonding, life usually re- 
quires reversible bonding. In both cases the strength of bonding is, therefore, a 
quantity of considerable interest. The most common measure of the strength of 
bonding is the work of adhesion A7, i.e. the energy per unit area needed to create 
two new interfaces when separating the two adhering objects. Many of the experi- 
mental methods developed to measure the work of adhesion actually rely on some 
combination of adhesive and elastic forces. For example, in the JKR test (Johnson 
et ai, 1971) the deformation of a soft sphere is used to infer the value of A7. Simi- 
larly, in the 'peel test', the work of adhesion is determined from the force required 
to peel a bonded film from a rigid substrate (Obreimoff, 1930; Kendall, 1975). 

As well as providing convenient techniques for measuring the work of adhesion, 
the interaction of elastic and surface forces gives rise to a rich variety of physical 
phenomena. In the peeling geometry, for example, fingering instabilies have been 
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observed for both viscous (McEwan & Taylor, 1966) and elastic (Ghatak et al., 
2000; Adda-Bedia & Mahadevan, 2006) substrates. Collective effects appear in the 
capillary- induced adhesion of an ensemble of flexible strips or rods (Bico et al., 
2004; Kim & Mahadevan, 2006; Py et al., 2007). These phenomena belong to a 
broad class of problems involving both elasticity and capillarity, which is referred 
to as 'elasto-capillarity'. 

Here, we consider a thin sheet adhering to a stiff substrate. In this situation, 
debonding can either be initiated from the periphery (as when peeling an adhesive 
strip) or from a point beneath the sheet. Both of these methods have commonly 
been used as a way of measuring the adhesive energy of polymer coatings on rigid 
substrates. In the latter situation, a shaft penetrating the substrate is used to 
push the sheet from below. This causes the formation of an internal delamination 
blister and so is commonly referred to in the literature as the blister test (Dannen- 
berg, 1961; Briscoe & Panesar, 1991). While this is a useful test for characterizing 
the properties of adhesives, it is distinct from the delamination observed in many 
manufacturing processes, which is often driven by pre-existing stresses within the 
substrate (Gioia & Ortiz, 1997). 

The experiment presented in this article is similar to the blister test except that 
the adhesion is mediated by the surface tension of a liquid — hence this experiment 
may be thought of as the 'liquid blister test'. More precisely, an elastic disk is bound 
to a substrate by a very thin liquid layer and quasi-statically loaded from below by 
means of a central indentor. This experimental setting allows us to go beyond the 
standard blister test and to unravel previously unreported regimes where the blister 
is not axisymmetric. We find an instability driven by orthoradial compression in an 
annulus, but the selection of the wavelength is different from that found in other 
systems (Mora & Boudaoud, 2006; Huang et at, 2007). The quasi-static nature 
of the experiment also distinguishes this instability from the fingering observed 
when a fluid displaces a more viscous liquid (Saffman & Taylor, 1958). At very 
large displacements, we find that the blister takes a triangular shape. The adhered 
part of the sheet then forms a conical shape reminiscent of the developable cones 
observed in crumpled paper (Ben Amar & Pomeau, 1997; Cerda et al., 1999). 

The paper is organized as follows. In §2, we detail the experimental setup and 
characterize experimentally the main regimes. These regimes are identified by the 
shape of the blister edge (circular, undulating or triangular) and depend on the 
displacement imposed by the indentor. In §3, we model the system as a thin elastic 
plate adhering to a stiff substrate with a constant work of adhesion. Using this 
model we study the properties of the blisters in each of the three regimes and 
compare the results to our experimental data. Finally, in §4, we summarize our 
findings and discuss some of their implications. 

2. The experiment 

(a) Experimental setup 

The system, see figure 1, consists of a thin elastic sheet adhering to a rigid glass 
plate (of dimensions 330 x300mm 2 ) by means of a thin liquid layer. The experiments 
reported here were performed with ethanol (surface tension 7 = 21mNm~ 1 ), but 
we also used silicon oils with no noticeable effect on the results. Ethanol was used 
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Figure 1. Experimental setup. Side (a) and plan (b) views, (a) A thin plastic sheet (see 
table 1 for properties) of radius R* adheres to a glass plate (330 x 300mm 2 ). The adhesion 
is mediated by a thin layer of ethanol. A hole of 16mm diameter is drilled in the middle 
of the plate for an indentor (cylinder of 6.3mm in diameter, capped by a hemisphere) 
that can produce displacements of up to 2cm. (b) When moving the indentor upwards, a 
delamination blister is nucleated and a crack front propagates (dashed curve). 



Material 


E (GPa) 


V 


h (/xm) 


B (mJ) 


2R* (mm) 


LDPE 


0.2 


0.4 


30; 100 


0.50 


120 


175; 275 


PP (i) 


2.1 


0.4 


40 


13 


120 


175; 275 


PP (ii) 


1.2 


0.4 


100 


120 


120 


175; 275 


PA 


2.7 


0.4 


200 


2100 


175 



Table 1. Properties of the sheets used: Young's modulus E, Poisson ratio v , thickness h, 
bending stiffness B and disk diameter 2R* . LDPE stands for Low Density PolyEthylene, 
PP for PolyPropylene, PA for PolyAmide-Nylon6 . The values of E were averaged over two 
orthogonal directions as these materials are anisotropic (anisotropy in the range 10-20%). 



for its ease of cleaning (compared to silicon oils) and the insensitivity of its surface 
tension to contaminants (compared to water). The plate had a hole of diameter 
16mm drilled in its centre, which allows for the equilibration of air pressure between 
the two sides of the sheet. We use a micrometric screw positioned beneath the 
hole as an indentor. Its cylindrical end is capped with a hemisphere and has a 
diameter of 6.3mm. The screw allows a vertical positioning of the cap to within 
1/xm over a 2cm range. The hemispherical cap of the indentor pushes the sheet 
causing debonding of the sheet and nucleation of the blister. As in the delamination 
literature, we shall use the term 'crack front' to designate the frontier between the 
bonded and the debonded parts of the sheet. Upon debonding, the sheet and the 
glass substrate remain wetted by the liquid, as shown schematically in the inset 
of figure la. The work of adhesion therefore corresponds to the energy required to 
create two air/ethanol interfaces of surface tension 7, and we have A7 = 27. 

Typically the indentor is moved at a vertical speed of lOO^ms -1 . Using typical 
speeds in the region of 5/xms _1 we see no observable differences. Likewise, with 
silicon oils of viscosities different to that of ethanol we do not see any observable 
difference. We therefore conclude that our system is effectively quasi-static and 
restrict the present study to equilibrium states. 

We used circular sheets made of three types of plastic material: Low Den- 
sity PolyEthylene (LDPE), PolyPropylene (PP) and Poly Amide (PA), supplied 
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Figure 2. The different regimes for the blister shape (sheets of diameter 2R* — 180mm). 
(a) A roughly circular crack front at small indentor height (d = 3mm, LDPE sheet, 
thickness h = 30/mi). (b) At moderate indentor height (d = 1.5 — 6mm), the crack 
front becomes unstable and oscillates with a centrimetric wavelength (PP, h — 100/im). 
(c) The undulation might evolve into a star shape (d — 10mm, LDPE sheet, thickness 
h = 30/im) . (d) The crack front becomes triangular after the tip of one 'finger' has reached 
the edge (d = 3mm, PP, h = 100/^m). (e) Sometimes the crack front becomes rectangular 
(d = 3mm, PP, h — 100/im), instead of triangular, (f) Schematics of the blister as seen 
from above. 

by Goodfellow. Their mechanical and geometrical properties are given in table 1. 
Typically, the Young's modulus E is of the order of lGPa, the Poisson ratio v « 0.4 
and the thickness h is in the range 30-200/^m, resulting in a bending stiffness 

Eh 3 , . 

^12(1^)' ( 2 -D 

in the range 0.5-2000^J. The diameter, 2R*, of the sheets was in the range 12- 
27.5cm. The initial volume of ethanol used as adhesive was in the range 20-1000/d 
depending on the radius of the disk. This corresponds to a liquid layer thickness in 
the range 1-50/xm. We found that most results were insensitive to the volume of 
ethanol, except in an intermediate regime that will be discussed below. We therefore 
did not systematically measure the quantity of ethanol used. Care was taken when 
laying the sheet on the substrate, in order to avoid pre-stress and to spread a 
uniform liquid layer. 

To aid visualization, the underside of the glass plate was painted white. The 
blister was lit uniformly using a halogen lamp. Top views of blisters were then 
taken using a digital camera, with a mirror inclined at 45° with respect to the glass 
plate. A home-made image processing method was used to detect the position of 
the crack front. 

Figure 2 shows images of the blisters as seen from the side. Three main regimes 
are observed and are shown schematically in figure 2f. At small indentor height d, 
the blister is roughly circular. At moderate d, the crack front oscillates. When the 
thickness of the ethanol layer is small, the undulations develop into a star shape at 
higher d. This undulating regime is the only one which is sensitive to the volume of 
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Figure 3. (a) Mean radius of a closed blister R vs. indentor height d. The symbols (see 
key) correspond to material type and sheet thickness as detailed in Table 1; all sheets had 
a radius R* — 175mm. (b) Top view and side schematic of a circular blister. 

the adhesive liquid. When d is further increased, one of the fingers of the debonded 
region reaches the edge of the sheet. The crack front then becomes an open curve, 
which is generally triangular, but sometimes rectangular. We focus in the following 
on the more generic regimes, i.e. circular, undulating and triangular blisters. 

(b) The mean radius of a closed blister 

The blister is approximately circular at small indentor height d and undulates for 
larger values of d. We measured the mean radius R of the blister as a function of d in 
both of these regimes. We found R(d) to be linear, with an offset of order 10/im for 
d = 0. This offset is due to a small indeterminacy in measuring the displacement 
when the indentor touches the sheet. Figure 3 therefore shows the curves R(d) 
shifted to ensure that they pass through the origin. The linear dependence R(d) also 
holds approximately when the blister crack front is undulating, and it is insensitive 
to the radius of the sheet R* (data not shown). Thus, the slope R'(d) appears to 
depend only on the thickness of the sheet, on its mechanical properties and possibly 
on the properties of the liquid. This dependancy will be clarified by the theoretical 
analysis. 

(c) The instability 

It appears that the blister loses axisymmctry for sufficiently large displacements 
d > d c . The threshold height d c at which the circular crack front starts to undulate 
is highly sensitive to both the thickness of the liquid layer and the preparation of 
the sheet on the substrate. d c tends to decrease for thick fluid layers or when the 
sheet is not perfectly flat initially. We were not able to determine this threshold 
d c in a reproducible manner. However, we were able to reproducibly measure the 
wavelengths of the undulations for the thinnest sheets. Figure 4 shows these ex- 
perimental results. We note that the wavelength of instability increases with sheet 
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Figure 4. The instability: undulations of the blister edge, (a) Top view of a portion of 
the crack front (LDPE, h = 30/Ltm, 2R* =175mm, d = 8mm). (b) Wavelength A of the 
undulations as a function of the indentor height d. In these experiments, we used 30/im 
thick LDPE sheets (diamonds) and 40/im thick PP (circles). For each material, we used 
two sheet diameters 27?* : 175mm (filled symbols) and 275mm (empty symbols). 

stiffness but seems rather insensitive to the radius of the sheet. Furthermore, the 
wavelength does not appear to be sensitive to the rate of loading and the undula- 
tions persist once the indentor stops moving — in contrast with the instability of 
Saffman & Taylor (1958). 



As the undulating crack front approaches the periphery of the sheet, a few un- 
dulations grow larger with a finger-like shape, as in figure 2c. With further increases 
in the indentor height d, one of the fingers reaches the edge first. This causes the 
other fingers to retract, and the crack front becomes an open curve. If the ethanol 
layer is not too thin (>20/im), the edge of the blister relaxes to a triangular shape 
rounded at one tip (figure 5a). This scenario is made possible by the visible sliding 
of the sheet along the substrate in the adhering region. Upon further increasing 
d, the angle at the tip of the blister at first increases but the tip disappears at 
sufficiently large values of d. The crack front then consists of two non-intersecting 
straight lines running between the sheet edges. We do not consider this regime here 
because of its ^reproducibility and instead focus on the range of indentor heights 
for which the rounded tip (B in figure 5a) is within a few mm of the edge of the 
sheet (C). 

We measured the angle (3 ~ 4a, defined in figure 5a, as a function of the inden- 
tor displacement d. This measurement is well-defined and insensitive to the curved 
shape of the tip. However, we plot the effective angle at the tip, 2a, as a function 
of tip displacement in figure 5b. In contrast to the earlier results for circular and 
undulating blisters, the radius of the sheet R* is an important parameter for tri- 
angular blisters. As a general trend, the opening angle increases when the bending 
stiffness B increases or when the radius R* decreases. We also measured the radius 
of curvature, r c , of the crack front at the tip as we shall use it in the theoretical 
section. Typical values of r c are in the range l-4cm. 



(d) The triangular blister 
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Figure 5. The triangular blister, (a) Top: plan view defining the angles a and /3. Bottom: 
schematic of a horizontal cross section at the tip of the blister, showing the meniscus (A), 
tip of the blister (B) and edge of the sheet (C). The distance AB is in the millimetric 
range, (b) Angle of the tip of the triangle 2a as a function of the indentor height d. The 
sheets were made of the three plastics (PP, LDPE and PA) with thickness in the range 
30 — 200/^m. Three disk diameters were chosen: 120mm (s), 175mm (m), 275mm (1). 

3. Theory and comparison with experiments 

The theoretical approach adopted here is to first parametrize the static blister shape 
in a given regime by a single parameter. For example, we use the blister radius, R, 
as the relevant parameter for circular blisters and the apex angle, 2a, for triangular 
blisters. This allows us to calculate the deformation field in the elastic sheet using 
thin plate theory (Mansfield, 1989). Once the deformation field is known, the elastic 
energy of deformation can then be calculated as a function of this parameter. Adding 
the surface energy AjA, where A is the area of the blister, we obtain the total energy 
of deformation £. The preferred value of the shape parameter, e.g. R or 2a, may 
then be determined by minimizing the energy £ with respect to variations in this 
parameter. 



(a) The circular blister 

(i) Scaling analysis 

We begin by considering a scaling approach for determining the properties of 
axisymmetric blisters. This approach will be verified in the more detailed analy- 
ses that follow, which will also yield the various prefactors in these scalings. An 
axisymmetric blister of radius R with a vertical displacement d at its centre is 
subject to a stretching strain e ~ d 2 /R 2 from elementary geometry. The sheet 
has a stretching stiffness (or spring constant) Eh and so the stretching energy 
£ e ~ Ehe 2 R 2 ~ Ehd 4 /R 2 . The bending energy Sb is proportional to the curvature 
(~ d/R 2 ) squared, so £ b ~ B(d/R 2 ) 2 R 2 ~ Bd 2 /R 2 . Finally, the surface energy 
required to create the blister is simply the energy required per unit area, A7, mul- 
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Figure 6. Setup and notation for the analysis of an axisymmetric blister. 



tiplicd by the blister area and so we have £ s 
system then takes the form 



A7_R 2 f. The total energy of the 



£ 



Bd 2 Ehd A 



R 2 



R 2 



A 7 i? 2 . 



(3.1) 



In (3.1), the blister radius is the only unknown and so we minimize this ex- 
pression with respect to R to find the blister radius B that is energetically the 
most favourable. There arc two limits that are of particular interest depending on 
whether d <C h or d ^> h. When d <C h, £& ^> £ e and so we find 



R 



a 7 ; 



1/4 



Conversely, when rf> h, £ e ^> £ b and we find 

1/4 



R 



a 7 ; 



(3.2) 



(3.3) 



Our experiments are conducted with d h and so we expect to observe the 
scaling (3.3). The argument above holds as long as the stretching energy does not 
vanish. In particular, it applies to axisymmetric configurations which necessarily in- 
duce stretching: an axisymmetric sheet takes the form of a cone of revolution, which 
is impossible for an initially flat sheet unless it is allowed to stretch or is cut. How- 
ever, as soon as the axisymmetry is broken, configurations with no stretching energy 
become possible (they correspond to developable configurations, see Ben Amar & 
Pomeau, 1997; Cerda et at, 1999; Cerda & Mahadevan, 2005, for example). For 
such situations, these scaling analyses reveal that stretching is enormously expen- 
sive (relative to bending) and so we expect the system to accommodate deformation 
via bending. We now move on to determine the prefactor in the scaling relationship 
(3.3) for an axisymmetric blister. 



(ii) General theory 

A theoretical idealization of an axisymmetric blister is shown in figure 6. In this 
geometry, the Foppl-von Karman equations governing the vertical displacement w 

f The meniscus near the contact region extends over a region on the order of the capillary 
length £ c = (7/pg) 1 / 2 where p is the liquid density and g the acceleration due to gravity. The 
gravitational potential energy of the liquid stored in the meniscus is then S g ~ pgt 2 c x Ri c ~ ^m c . 
Provided that R 3> £ c we find that E s ^> £ g . This is always the case in our experiments and so 
we neglect this effect in everything that follows. 
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of the sheet as a function of the radial co-ordinate r become (see Mansfield, 1989, 
for example) 

A f A V 1 A 1 1 a 

(3.4) 



B—\r— 
dr 1 dr 



and 



dr I dr 



1 d 

r dr 

1 d 

r dr 



,rw 



} - -Kw 



(3.5) 



Here the bending stiffness £? is given in (2.1), q is a load distribution and the Airy 
function is a potential for the stress: 

a rr = —r- and a ee = TT- ( 3 - 6 ) 
r dr dr A 

For a blister of radius R, we introduce dimensionless variables 

p = r/R, W = w/R, Q = qR/Eh, $ = 4>/EhR 2 . (3.7) 

Within the blistered region (p < 1) we assume a point force acts at the origin so 
that Q(p) = Q5(p)/2irp. 

Substituting the dimensionless variables in (3.7) into (3.4) and (3.5) and inte- 
grating once with respect to p we find that (Jensen, 1991) 



dp 



and 



where 



Id, f \ 



B 



2tt 



= -±W' 2 

2 ' 



1 



EhR 2 12(1 



i? 2 ' 



(3.8) 



(3.9) 



(3.10) 



For the values of h and R relevant in our experiments, e <C 1. We therefore ex- 
pect that the effects of bending may be neglected in comparison with stretching. 
Furthermore, this assumption should be valid throughout the sheet except in small 
boundary layers near p = 0, 1. 

As already discussed, bending dominates stretching for small vertical displace- 
ments and there is a transition between bending and stretching behaviours at in- 
termediate displacements, as our earlier scaling analysis shows. Cotterell & Chen 
(1997) consider a related problem but their analysis is valid only for d/h < 10. The 
transition between bending and stretching has also been studied by Wan (1999), al- 
beit under the assumption that the radial and tangential stresses in the membrane 
are equal and constant. In Appendix B we consider the two-dimensional version of 
this problem. We show there that the scaling (3.3), with the appropriate prefactor, 
is correct to within 10% once d > Ah. The assumption that d » h is true for all 
but the very early stages of our axisymmetric experiment and so we shall consider 
here the 'membrane limit' e = 0. 

Setting e = in (3.8) we find 



Q_ 

27T' 



(3.11) 
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which is to be solved along with (3.9) and the boundary conditions 

W(0)=d/R = W o , W(1) = 0, $'(0) = 0. (3.12) 

(We note that we may arbitrarily set $(0) to zero because only $' enters our equa- 
tions.) A fourth boundary condition is required to close the system. This boundary 
condition should be on the horizontal displacement, 

U r = -j^ivee ~ Wrr), (3.13) 

at the edge of the blister where o~ rr and agg are the radial and circumferential 
tensions in the sheet and are given by (3.6). However, there are two plausible 
conditions on u r (R). If the sheet does not slide outside the blister, then u r {R) = 0. If 
sliding is permitted then we also need to solve the equilibrium equations outside, and 
we require only continuity of displacements at the blister edge, u r (R + ) = u r (R~). 
In practice, we might expect some combination of these conditions to be realized. 
We shall therefore consider both of these possibilities in the following. 



(iii) No sliding 

If the horizontal displacement at the blister edge is constrained to be zero, the 
constitutive equation (3.13) and the definition of the Airy stress function (3.6) yield 
the fourth boundary condition 

$"(1) - = 0. (3.14) 

We look for solutions of (3.9) and (3.11) of the form W = W (l - p a ), $ = bp c , 
and find that 

w = w (i- P 2 / 3 ), * = %w* P v 3 (315) 

This solution is well known and is given, in slightly modified notation, in Mans- 
field (1989). However, this solution docs not satisfy the boundary condition in (3.14) 
unless v — 1/3. For values of v « 1/3, a perturbation of (3.15) allows us to satisfy 
this boundary condition, as shown in Appendix A. However, the perturbative so- 
lution only changes the final result of the analysis (the predicted blister radius R) 
by less than 3% for v — 0.4, which is a typical value in our experiments. We shall 
therefore use only the solution in (3.15) in the main text. 

The deformation field in (3.15) gives rise to dimensional stresses (see Eq. 3.6), 

arr = Eh ^p-V\ aee = EhM-p-^ (3.16) 



in the sheet. The stretching energy within the blister is then given by 

\ J o- a pe a p dA = -^j- J r{a 2 rr - 2vo rr o ee + a 2 e ) dr 

tt(5 - Zv) Ehd A 
48 #2"' 



(3.17) 



where the components of the stress tensor, a a p, are related to the components of 
the strain tensor, t a p, by the linear elastic constitutive relation (Mansfield, 1989). 
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Combining the stretching energy in the sheet with the change in surface energy 
associated with a blister of radius R, we find that the total energy of the system is 

£ = ^E^_ + ^ R2 (318) 
Choosing R such that dE/dR = then gives that 

R =h{(l--)^} 1/4 d- (3-19) 



(iv) Sliding 

If sliding is allowed beyond the blister, we only require the continuity of hori- 
zontal displacement at p = 1, 

$"(1+) - = $"(1-) - (3.20) 

Within the blister, it is convenient to introduce 

^)=p^{p), V = p 2 . (3.21) 

Along with (3.8), this substitution simplifies (3.9) to 

^--^V 2 (3 22) 

which integrates to 



dip Q {Aip + l^ 1/2 



drj Air \ ip 



(3.23) 



We then consider the region outside the blister (p > 1) where, by assumption, 
we have that W = 0. However, there is a non-trivial stress function <J> in this part 
of the sheet. The most general solution of (3.9) for $ is $' = ap + b/p. Requiring 
that <7 rr (oo) = gives a = 0. The constant b = by the continuity of a rr at 
p = 1, so that 

<&' = </>(!)//>• (3.24) 

The continuity of displacement (3.20) then gives = 0, which in turn determines 
A = -1/V>(1) in (3.23). Integrating (3.23) analytically, subject to $'(0) = 0, we 
find that 

l^-y^l^+tan^^L, (3.25) 

where tp = (We also note that the requirement that V(l) = 1 yields 

Q = 47nA(l) 3/2 .) This equation gives ip implicitly in terms of rj. We may also 
integrate (3.8) to obtain W as a function of tp, rather than r\. We find 



W = W 



i-isin- 1 (2^-l) , V(l) = WS/n 2 . (3.26) 
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Figure 7. The axisymmetric blister. Experimental results for the non-dimensional mean 
radius of the blister R/R* as a function of non-dimensional rescaled indentor height 
d/[R* (Eh/ A'y) 1 / 4 ']. Here E is Young's modulus, A7 is twice the surface tension of ethanol, 
h is the thickness of the sheet and R* its radius. The theoretical line predicted by model 
1 (no sliding outside the blister) has a slope of 0.53 (solid line), while model 2 (allowing 
the outside to slide) gives a slope of 0.47 (dashed line). 

The stretching energy for the total deformation of the sheet may be calculated 
analytically to be 

1 Ehd 4 

^=2^, (3-27) 

independently of the Poisson ratio v. Adding this to the surface energy and mini- 
mizing the total energy with respect to variations in R, as in the previous section, 
we find that 



l_Eh\ lf \ 
2tt 2 A 7 ) 



R= 



(3.28) 



(v) Experimental results 

We first note that equations (3.3), (3.19) and (3.28) all reproduce the linear 
dependence of the blister radius R on the indentor displacement d observed in 
experiments. These models also predict that the constant of proportionality is itself 
proportional to s = (Eh/ A7) 1 / 4 . To go further, we rescale the experimental data of 
figure 3 using the prefactor s. This rescaling allows a collapse of the data as shown 
in figure 7. The slope of the experimental data is 0.5 ± 0.1, which lies between the 
slopes (5/3 - vf/ i /2 « 0.53 and (2tt 2 )- 1 /4 ~ .47 given by model 1 (no sliding 
outside the blister) and model 2 (sliding allowed), respectively. The vertical scatter 
in the rescaled data are most likely due to the relatively large errors in the value 
of Young's modulus (10-20%) caused by the inhomogeneity and anisotropy of the 
sheet. Although we could not visualize sliding in circular blisters, we believe that it 
is more likely to occur when the contact is lubricated by a relatively thick ethanol 
layer; solid friction would prevent sliding if the ethanol layer is very thin. This 
observation has implications for the instability of the circular blister as discussed 
in the next section. 
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The scaling in (3.19) and (3.28) has been observed previously in the conventional 
blister test (Briscoe & Panesar, 1991) with a theoretical analysis presented by Wan 
& Mai (1995). However, we note that the coefficient there is slightly different since 
the analysis of Wan & Mai (1995) explicitly assumes that the beam deflection takes 
the form of a cone W = W (l — p). Here we have determined both W and from 
the fully coupled problem. 



(b) Instability 

(i) Theory 

In the absence of sliding, we are unable to propose any instability mechanism. 
However, when sliding is allowed, we find that the azimuthal stress, as resulting 
from (3.24), 

a e0 = -Eh/ir 2 {d/rf < 0, (3.29) 

is compressive outside the blister. (Note that there is always azimuthal traction 
inside, see e.g. Eq. 3.16.) This compressive stress could lead to a buckling instability 
that we analyse here. 

We suppose that this instability gives rise to a vertical displacement of size H 
with wavelength A, which we assume decays over some radius R. The bending energy 
of these fluctuations should be of size BH 2 R 2 /X 4 . A compressive energy within the 
sheet is released by going out of the plane. The perturbation to the stretching 
energy has typical size ~\agg\H 2 R 2 /\ 2 . Finally the undulation of the crack front 
incurs an energy penalty. The amplitude of the contact line oscillations is H/Wq 
where Wo is the slope of the blister, and we find an energy penalty Aj(H/Wo) 2 
coming from the perturbation to the area of the blister. Adding these three energy 
terms, we find that 

,^(£_M + _^). (3 ,0, 

We see that if the term in brackets is positive then increasing the perturbation 
amplitude H increases the total energy and so is discouraged: the system is stable. 
Conversely if this term is negative, energy may be released by increasing the per- 
turbation amplitude: the system is unstable. For the system to remain stable we 
therefore require that 

M ~£ + #f§ s « A >- < 3 - 31 > 

The minimum value of the function /(A) is attained with 

,.^WV /4 -(^V /4 m'74V" (3.32) 



The equation (3.31) then requires that \agg\ < min/ for stability, i.e 



BA-f 



1/2 
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^ l 'J^Li 

100 200 300 
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Figure 8. The instability. Non-dimensional wavelength X/[h(Eh/ A7) 1 / 4 ] as a function of 
the non-dimensional indentor height d/h. The data and symbols of figure 4 are replotted 
here. The data collapse on a master curve which can be fitted by a power law (inset and 
solid curve) with an exponent of 0.32 ± 0.05. The dashed curve shows a best fit to the 
predicted square root dependence. 



Using the expression for age in (3.29) evaluated at r = R (where the compression 
is largest) we may write the condition of stability as 

d/R<h/R. (3.34) 



This scenario involves fluctuations of the sheet on both sides of the crack front. 
The compressive stress ago outside the blister is comparable to the stress inside as 
long as r ~ R (3.29), suggesting that the driving mechanism for the instability acts 
on a scale i? ~ R. The condition (3.34) then implies that the system is unstable as 
soon as d > h. 



(ii) Experimental results 

As seen above, the system is unstable whenever sliding occurs. This accounts 
for the sensitivity of the experimental threshold to the preparation of the sheet and 
the thickness of the liquid layer: sliding is facilitated by the lubrication of a thicker 
layer. 

Following (3.32), we rescale the experimental data for the wavelength by plotting 
\/[h(Eh/A^y/ 4 ] as a function of d/h. Figure 8 shows a satisfactory collapse of the 
data onto a master curve. The power-law dependence is found to be 0.32 ± 0.05, 
which is weaker than the predicted 1/2 power. This might be attributed to the 
approximation R = R. In fact, R can be smaller than R when the stress field feels 
the outer boundary of the sheet — an effect neglected when we determined the 
stress function (3.24). Note that other mechanisms of wavelength selection (such 
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as those studied in Mora & Boudaoud, 2006; Huang et al, 2007) would lead to no 
dependence on the Young's modulus and so cannot account for our results. 

(c) Triangular blisters 

(i) Theory 

Our earlier analysis of a circular blister was based on neglecting the energy due 
to bending. This simplification was valid for two reasons: firstly the energetic cost 
of bending is much smaller than that of stretching and secondly an axisymmetric 
configuration necessarily induces stretching. When the crack front becomes an open 
curve, the system has more freedom and thus can avoid expensive stretching. Our 
analysis for triangular blisters assumes that the sheet does not undergo in-plane 
stretching and so has a developable conical shape (see also Ben Amar & Pomeau, 
1997; Cerda et al., 1999; Cerda & Mahadevan, 2005). In our system, the adhered 
part of the sheet is free to slide eliminating the possibility of a 'hoop stress', which 
is an additional complexity in previous analyses (see Cerda & Mahadevan, 2005, 
for example). For simplicity we derive here the governing equation for the deformed 
shape of the cone, though it may also be recovered as the small deformation limit 
of equation (4.17) of Cerda & Mahadevan (2005) in the absence of a constraint, 
i.e. a = 1. 

A theoretical idealization of a triangular blister is shown in figure 9a. Here 
we treat the angle at the tip of the triangle, 2a, as the unknown parameter to be 
determined. We use a cylindrical polar co-ordinate system, (r, 9), centred on the tip 
of the isoceles triangle. (Note that r is dimensional, since the non-dimensionalization 
used in earlier sections is now redundant.) The only developable shape compatible 
with a triangular crack front is a cone, so that 

w(r,6) = d^f(6) (3.35) 

for some unknown function f{9). Here, we shall only consider 9 > — the symmetry 
of the sheet shows that in fact / depends only on \9\. 
The bending energy of the sheet is then given by 

r-a r2R* 

S b ss B / d6 / dr r(V 2 w) 2 

= f^/V + />^, 

where r c is some cut-off length that we assume is unimportant in the following. 

Requiring the bending energy in (3.36) to be stationary with respect to varia- 
tions in /, we find that 

f(iv) + 2 /(«) + / = 0. (3.37) 

Equation (3.37) has solutions of the form f(9) = (A6 + B) sin6> + (C9 + D) cos9 
where A, B, C, D are constants to be determined from the boundary conditions, 
which correspond to continuity of displacement and slope: 

/(0) = 1, /'(0) = 0, f(a) = 0, f(a) = 0. (3.38) 
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Figure 9. The triangular blister, (a) Theoretical idealization of the triangular blisters 
observed at large displacements, (b) Experimentally measured angle at the tip of the 
triangle, 2a, as a function of 8 — Bd 2 /[AyR* 4 ] log 2R* /r c where B is the bending stiffness, 
A7 = 27 is twice the surface tension of the liquid, R* is the radius of the sheet and r c is 
a cut-off which corresponds to the radius of curvature of the tip. We measured /3 — 4a for 
three types of plastic material (PP, LDPE, PA) with thicknesses ranging between 30^im for 
LDPE and 200/^m for PA. Three sizes of disk were used: 120mm (s), 175mm (m), 275mm 
(1). The solid curve shows the theoretical prediction obtained by minimizing (3.42). 

Thus the cone profile is given by 

(a + sin a cos a) (sin \0\ - \0\ cos |0|) - flsinflsin 2 a 

}{O)=cos0-\ ;-= . (3.39) 

ar — sin a 



We shall adopt the same approach as earlier writing down the total energy of the 
system and minimizing this with respect to variations in the angle a. The bending 
energy is given by 

. Brf 2 2a + sin2a 1 2R* 

£ b~-^2 — • 2 l0 § — ■ ( 3 " 4 °) 

R* or - sin a r c 

Using elementary geometry, the change in surface energy of this configuration rel- 
ative to the flat state is 



£ s = A7 J R* i (2a + sin2a), 



(3.41) 



where R* is the radius of the sheet. 

Adding the two energies in (3.40) and (3.41), we find that the angle a should 
be chosen to minimize 



g(a) = 2a + sin 2a + 6 



2a + sin 2a 
a z — sin a 



(3.42) 



where 



Bd 2 , 2R* 
4 log 



A 7 i? 



(3.43) 



For minima we require g'(a) = 0, which can easily be solved numerically, yield- 
ing the solid curve shown in figure 9b. For a <C 1, the leading order behaviour of 
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g(a) may be used to show that the total included angle at the apex of the triangle 
(measured in degrees) is 



We note that this scaling is consistent with the scaling analysis of §3 (a)i: replacing 
R with aR* in (3.2) we recover (3.44). Furthermore, we find that (3.44) is correct 
to within 5% provided that 5 < 0.01. 

(ii) Experimental results 

In figure 9, we plot the angle at the tip of the triangle as a function of the 
non-dimensional variable 5 defined by (3.43). We took the cut-off length r c to be 
the radius of curvature of the crack front at the tip of the triangle. There is good 
agreement with the theoretical curve, showing that this regime is ruled by a balance 
between bending and capillary forces. 



In this article we have considered a thin sheet adhering to a stiff substrate by 
means of the surface tension of a liquid. We have investigated, experimentally and 
theoretically, the debonding of the sheet under an imposed central displacement. 
This is the generic debonding mode whenever the sheet is not peeled from an edge. 
We uncovered three main regimes dependant on the central displacement. For small 
vertical displacements, the debonded area is axisymmetric. This regime was known 
previously in the context of the conventional blister test (Dannenberg, 1961; Briscoe 
& Panesar, 1991) where there is no liquid. The presence of a liquid in our system 
allows sliding and so we derived a new analytical solution to account for this. In 
the second regime, the crack front undulates, and we determined a characteristic 
wavelength for the undulations. Finally, the crack front becomes triangular, and 
the shape of the sheet is approximately a developable cone. Overall, our theory is 
in good agreement with experiments suggesting that the present framework might 
serve as a reference for the adhesion of clastic sheets. 

Our study may also be relevant to small scale biological and technological sys- 
tems such as the hair of gecko (Autumn et al, 2002) or microelectromechanical 
devices (Xia et at, 1999; Mastrangclo & Hsu, 1993): capillary and van der Waals 
forces become increasingly important at small scales. The delamination of thin films 
in general is slightly different than that studied here. In general settings, debonding 
is driven by in-plane pre-stress and so there is a compressive azimuthal stress in the 
debonded part also. However, we expect our result for the wavelength in the un- 
dulating regime to hold for the stability of delamination blisters, which apparently 
has not been investigated previously. Despite this expectation, much more work is 
required to fully understand the complex patterns generated by the competition 
between elastic and adhesive forces. 

We are grateful to Mokhtar Adda-Bedia for fruitful discussions. D.V. is partially supported 
by the Royal Commission for the Exhibition of 1851. 




(3.44) 



4. Conclusion 
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Appendix A. Membrane solution with v ^ 1/3 

In this Appendix we consider again the membrane solution for the blister given in 
§3a. Recall that (3.15) was used as the solution of the membrane equations but 
does not satisfy the no-slip boundary condition at the blister edge, (3.14), unless 
v = 1/3. Here, we look for a small perturbation to the approximate solution (3.15) 
of the form 



W 



W (l-p 2/3 )+SW(p), $«^^ V /3 + ^(p). (Al) 



Substituting these forms into the membrane equations (3.9) and (3.11) and assum- 
ing that SW, i$ C 1 we find that 



W « Wo (l - P 2/3 + Av[p 2 - p 2 / 3 ]) , $ „ ±W 2 p 4 / 3 (l + Ai/[2 • 



§P 4/3 ] 



(A 2) 
where 

Typically in our experiments v — 0.4, which corresponds to Av w 0.02 justifying 
our assumption of a small perturbation. 

Using the value of $ in (A 2) to calculate the stretching energy in the sheet we 
find, to leading order in Av, that 



•4 



£e^(f--)^(l + 6A,). (A4) 

Performing the same total energy minimization as earlier we find that the radius 
of the blister is 

*«i{(§-")|£} d[l + lAv + 0(Av 2 )] (A5) 

Using v = 0.4 we find that the correction to the blister radius provided by (A 5) of 
the earlier expression (3.19) is < 3%. 



Appendix B. Analysis of a Two-Dimensional Blister 

In this Appendix, we consider the two-dimensional version of the axisymmctric 
blister test treated in the paper. In particular, we shall determine the effect of a 
non-negligible bending stiffness on the radius of the blister for a given indentor 
height and quantify more precisely when it is reasonable to use the membrane 
approximation. We use a two-dimensional geometry since it is then possible to 
make considerable analytical progress. 

For a two-dimensional geometry we find that the Airy stress function is given 
by (j) — a xx y 2 /2 for some tension o xx . The first Foppl-von Karman equation then 
takes the form 

d 4 w d 2 w 
B d^ = a **dx*- (B1) 
Using the non-dimensionalizations defined in (3.7) of the paper and letting 

" = £ Eh (B2) 
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wc find that 

d±W _ d 2 W 

dX^-^d^' (B3) 
which is to be solved with the continuity of slope and displacement 

W(0)=W , W'(0) = 0, W(l) = W'(l) = 0. (B4) 

The solution of (B 3) with boundary conditions (B 4) is 

, , (1 - |X|)cosh±u + «~Hsinh/i(|X| - \) -sinh±u| 

W(X) = W - — ^—f- — l - PI ' 2 ' (B5) 

coshi,u — 2[i i sinh^/x 

Determining the value of a xx , and hence /j,, requires an additional condition. This 
condition is the relationship between the strain e xx and the tension a xx (Mansfield, 
1989) 



dU ,fdW\ 2 o xx (\-v 2 ) 

dX + 2 {dx) ~ Eh ' (B6) 



Integrating this relationship between X = and X = 1 and assuming that the 
sheet is fixed (£7(0) = U(l) = 0) at the boundaries shows that 



2a xx (l -v 2 ) 
Eh 

where 



f 1 W'(XfdX = Wlf 1 (p), (B7) 
Jo 



/^e 2 + cosh A .-3^sinh^ (Bg) 
2 (cosh i^ — 2\i 1 sinhi^i) 

(We note that when a two-dimensional sheet is free to slide, there will be no stretch- 
ing and only bending of the sheet. This is qualitatively different from the axisym- 
metric case where, even in the presence of sliding, geometry forces the sheet to 
stretch.) 

For our purposes, it is useful to eliminate a xx from (B 7) in favour of /i by using 
(B2). Introducing 

"■A"'?' (B9) 

we find that \i satisfies 

M 2 /u; = /i(M), (BIO) 

which has solution \i = /Lt c (w). (The advantage of using u> rather than Wq/c is 
that uj is independent of the blister radius i?, which is as yet unknown, lo measures 
the vertical displacement relative to the sheet thickness and so is a natural control 
parameter.) 

We continue by considering the limits of small and large u>, which, as we shall 
see, correspond to small and large values of /i, respectively. From the asymptotic 
behaviour of /i(/i) for small and large fi we find that 

lfc R "' / x , ^, ( B11 ) 

1, u> » 1. 
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Note that the two cases w < 1 and w»l give rise to only very slight variations in 
the asymptotic relations for /j, c . 

To determine the blister radius, R, we must calculate the deformation energies 
of the sheet. The bending energy is given by 



£t = b j\"{xf Ax = ^/ 2 (m); Mm) = § 



sinh /i — /i 



2 (cosh |/x — 2/i _1 sinh ^fij 

(B12) 

The stretching energy is given by 



£ e = o xx \ w\xf dx = EhReW^'hip), (B 13) 

Jo 

with fi(fi) as defined in (B8). We can therefore write the total (elastic) energy as 
£/{EhR) = eW 2 [/ 2 (/x c ) + /ic/i(Mc)] - 2e 2 (l - i/ 2 )wffM (B 14) 

where 

</(w) = / 2 [/i c (w)] + Mc(w) 2 /i [/ic(w)] • (B 15) 

and n c is given by the solution of (B 10) for a given value of ui. 

Finally we can now include the surface energy to determine the size of the 
blister. The surface energy gained is 2AjR so that the total energy of the system 
is 

£ = 2^R + 2{l-v 2 )EhRe 2 ujg{uj). (B 16) 

Recalling that e oc R~ 2 we may then minimize this expression by varying R. We 
find that the value of R at which the energy is minimized may be written 



V4 /„(, ,\ \ 1/4 

(B17) 



R _ ( 3Eh y* f g(u) 

The function {g{ui) /w} 1 ^ 4 , which determines R/d via (B 17), is plotted in figure 
10. We note that as u — > exo, g{uj)/u — > 1 so that 



R ( 3Eh 



d V 4A 7(!-^ 2 ) 



(B18) 



and we obtain the same scaling relationship as for the axisymmetric analysis. We 
observe that ui > 100 is sufficient for the relationship (B 18) to hold to within 10%. 
Furthermore, this corresponds to d > Ah — a condition that is easily satisfied in our 
experiments. We therefore expect that our neglect of bending effects in the analysis 
of §3(a) should be valid except in the earliest stages of the experiment. 

For completeness, we note that for w < 100 the prefactor is bigger and indeed 
that the scaling R ~ d no longer holds. For w « 1 wc find that g(u) w 12 so that 



R ( 3Eh ^ 1/4 



d \2Ajil-u 2 ) 



(B19) 
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Figure 10. The function {^(a;)/^} 1 ^ 4 plotted on logarithmic scales. From (B 17) we see 
that this function determines the evolution of the blister aspect ratio, R/d, as d changes 
with the material properties of the system fixed. 




Figure 11. Non dimensional blister shapes for different values of the parameter u>, defined 
in (B9). Shapes are shown for u> = 1 (dotted curve), u> — 10 2 , w = 10 3 and uj — oo (solid 
line) . 
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